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Abstract 



The formation process of a Bose-Einstein condensate in a trap is described 
using a master equation based on quantum kinetic theory, which can be well 
approximated by a description using only the condensate mode in interaction 
with a thermalized bath of noncondensate atoms. A rate equation of the form 
h = 2W~^{n) |(l — e{Mn-Ai}/fcT~j 77, + l| is derived, in which the difference be- 
tween the condensate chemical potential Hn and the bath chemical potential 
fj: gives the essential behavior. Solutions of this equation, in conjunction with 
the theoretical description of the process of evaporative cooling, give a char- 
acteristic latency period for condensate formation and appear to be consistent 
with the observed behavior of both rubidium and sodium condensate forma- 
tion. 
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Typeset using REVTeK 



The experiments on Bose-Einstein condensation of dilute atomic gases [|T]-^ have stimu- 
lated theoretical effort, which has however not produced any definitive result for the growth 
of the condensate from the vapor, although there have been significant theoretical contri- 
butions This paper will present a quantitative and experimentally testable description 
of the growth process, based on quantum kinetic theory P,p!0| , which can be simplified to a 
single first order differential equation for the number n of atoms in the condensate. 

Our formulation contains the following principal features. We use the Hamiltonian 

+ d^X J d3xV^(x)^^(x')M(x-x')^(x')^(x) 

+ J rf=^x\/T(x)^t(x)V'(x). (1) 

The potential function ^(x — x') is as usual not the true interatomic potential, but rather a 
short range potential — approximately of the form m5(x — x') — which reproduces the correct 
scattering length. 

We divide the condensate into two regions called the condensate band Rc, and the non- 
condensate band Rnc, as in Fig.l. We treat Rnc as being thermalized, representing the 
majority of the atoms as a heat bath which provides the source of atoms for condensate 
growth. The condensate band is the region of energy levels less than a value En, which in- 
cludes not only the ground state, in which the condensate forms, but also those levels which 



would be significantly affected by the presence of a condensate. |T3| In the noncondensate 
band, with energy levels greater than Eji, there is no significant such effect. 

The behavior in Rc is treated fully quantum-mechanically, and a description in terms of 
trap levels modified by the presence of a condensate is used. At any time there is a given 
number N of atoms in Rc, and the energy levels in such a situation can be described using 



the number-conserving Bogoliubov method devised by one of us [|T4[, so that the state of 
Rc is fully described by the total number of atoms in Rc, and the quantum state of 
the quasiparticles within Rc- In this formulation we can write the condensate band field 
operator in the form 



^c(x)=i?(eiv(x) + E 

I. m 



N 



(2) 



The quasiparticles, of energy e]}, are described by annihilation operators bm, while is 
the creation operator which takes the Rc system, for any N, from the ground state with 
atoms to the ground state with A^ + 1 atoms. The condensate wavefunction is ^Ar(x), and 
this satisfies the Gross-Pitaevskii equation 

- T^V'^TV + Vt^n + Nu\^Nf^N = fiN^N. (3) 

2m 

The amplitudes /m(x), (7m (x) are for creation and destruction of quasiparticles of energy e^, 
and are fully defined in |lTO|,p^, but will not play any significant part in this paper. 
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In this number- conserving Bogoliubov method, the atoms are conserved, while the quasi- 
particles are mixtures of phonon states, and these phonons relate to the process of transfer- 
ring an atom from an excited quantum state to the condensate level. Thus, the operators 
bm, &m do not change the total numbers of particles, while the operator B, which multiplies 
everything else in (0), reduces the total number of particles by 1. 

The process we wish to describe is as follows: 

(i): Some of the collisions in Rnc will transfer an atom to Rc, so that N —>■ N + 1, 
and there is of course the reverse process where a collision of a noncondensate band atom 
with one within the condensate band transfers an atom from the condensate band into 
the noncondensate band, so that — > — 1. (ii): We consider a situation in which 
there is initially no condensate — however the boundary between Rc and Rnc is fixed to 
be appropriate for the amount of condensate which is finally formed, (iii): By evaporative 
cooling, the chemical potential of the atoms in Rnc becomes nonnegative; this is permissible 
provided the chemical potential does not exceed the lowest energy Eji of Rnc- (iv): With 
a weak interaction potential u the Bogoliubov spectrum and wavefunctions are valid for all 
N, large and small, since for small N and u the results are not significantly different from 
perturbation theory. 

Using quantum kinetic theory [|TU| it is possible to derive a simple master equation for 



the density operator p which describes the state of the condensate. The main processes 
are caused by an atom scattering into or out of Rc, and this can occur in six ways; that 
is — s> ± 1 with one of; no change in the number of quasiparticles, the creation of a 
quasiparticle or the absorption of a quasiparticle. 

The six transition probabilities can now all be written in terms of the functions R^ as 

(4) 
(5) 
(6) 
(7) 
(8) 
(9) 

The functions R'^{y,Lj) are defined by 

R^iy.io) = I d'^ I dTA{T,u)F,F2{l + Fs)Wy{^,k) (10) 

^-(l/,^) = 7^^ Jd'^JdTA{T,uj){l + F,){l + F2)FsWy{K,k) (11) 



W+{N) 


= R'^iC.N, 


Pn/K) 




= R-{^N- 


i,Pn-i/Tt) 




= R~^{fm, 


(e- + MA) 


W--{N) 


= R'ifrn, 


(e^-i+/ijv-i)A) 




= R^{gm, 


(-6^+ MA) 




= R~{gm, 


(-e^^_i+/ijv~i)A) 



(27r)5^2 

where we will use the notation 



dV = d^Kid^K2d^K3d^k (12) 

A(r, LO) = 5(ACU123(X) - UJ)6{K, + K2 - Kg - k). (13) 



Here we use the notation 



3 



to represent the Wigner function corresponding to the wavefunction ?/(x). The function 
Fi = F(Kj,x) is the noncondensate atom density per of phase space, and 

and At^i23(x) = t^Ki(x) + t^KaW - UJKsi^) (15) 

with hujKi-y^) = + Vt(x). (16) 

We can write a stochastic master equation for the occupation probabihties p(A^, n), where 
n = {nm}, the set of all quasiparticle occupation numbers, in the form 



p{N, n) = 2NW+{N - l)p{N - 1, n) - 2(A^ + l)W+{N)p{N, n) 
+ 2{N + l)W-{N + l)p{N + 1, n) - 2NW-{N)p{N, n) 
+ X:{2n^M/++(iV - l)p{N - 1, n - e„) - 2{n„, + l)W^++(iV)p(iV, n)} 

m 

+ 5]{2(n^ + l)W--{N + l)p(iV + 1, n + e„) - 2n„W^„-(iV)p(iV, n)} 

m 

+ E{2(^m + l)W^+^(iV - l)p(iV - 1, n + e„) - 2n^l^+'(iV)p(iV, n)} 

m 

+ 5]{2n„W^-+(iV + l)p{N + 1, n - e„) - 2(n„ + l)H^-+(iV)p(iV, n)} (17) 

m 

Here = {. . . 0, 0, 1, 0, 0, . . .} has its only nonzero value at the position corresponding to 
the index m. 

We can interpret functions R^iy^ui) as forward and backward collision rates for those 
collisions which result in a particle entering (+) or leaving (— ) the condensate with an 
energy huo. The collision must also take place in a position where the condensate Wigner 
function is nonzero. The terms in (p!7[) representing transitions to the ground state of the 
condensate exhibit a stimulated increase in collision rate of approximately N , which can be 
a number up to 10'', but the transition probabilities defined in (|^-|^) are multiplied 

only by n^, which does not become large. Thus as an initial approximation we drop the 
quasiparticle terms, which are smaller by a factor of N , and are thus negligible for most of 
the condensation process. We note that the condensate wavefunction is in practice sharply 
peaked at x = by comparison with the phase space distribution function F(K, x) , and 
thus replace x wherever it occurs by 0, except in VFg^(x, k), whose integral gives the k space 
condensate probability density |^Ar(k)p. We finally get the simple master equation which 
consists of only the first line of ([l^); and in which the transition matrix elements take the 
simplified form 

W^{N) = j^-^ J dTA{T,^^{N)/n) 

xFiF2(l + F3)|f^(k)|2 (18) 
W-iN) = j^-^ J dTAiT,^^iN)/n) 

x(l + Fi)(l + F2)F3|eAr(k)|2 (19) 
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The evaluation of can be done in various degrees of approximation; we shall take here 
a thermal Rjyc with 



F{K, x) ^ ^e^^-^^+^Ti^^-'^'^/kT _ (20) 

from which one easily obtains (choosing Vt{0) = 0) 

W-^{N) = e^^'-''''^/''^W-{N), (21) 

which indicates that a steady state is reached at large when the chemical potential of the 
condensate almost equals that of the noncondensate. (Inclusion of the smaller terms in the 
master equation also shows that the temperatures of Rc and R^c are equal in equilibrium). 
Evaluation of can be done by taking the energy range of Rc to be negligible compared to 
that of -F(K, 0), and by taking the range of k to be small compared to that of K in -F(K, 0). 
We also approximate the Bose function ( pOD by its Boltzmann equivalent for most of the 
range of integration since the integrals can then be evaluated analytically — however this is a 
purely technical issue, which does not affect the essence of the results. Using u = Anah'^/m, 
where a is the s-wave scattering length, we get 

WHN) = 'j^^e"-'^-{p,,{!^)}. (22) 

(Here Ki{z) is a modified Bessel function. Notice also that the prefactor Am(akTy / nh'^ is 
essentially the elastic collison rate pav , where the quantities are evaluated at the critical 
point for condesation.) Under the assumption that the majority of the atoms are in the 
condensate, the major behavior of the master equation ([T7|) is given by the rate equation 
for the mean number of atoms in the condensate (written as n), 



n 



2W+{n) { (l - e^^"-^>/'^^) n + l}. (23) 



Since ( PT| ) has been used, this represents a situation in which a condensate (which may 
be initially unoccupied) is in contact with a bath of noncondensed atoms. If (^3]) is used 
alone, it describes a situation in which the thermal bath is not depleted as the condensate 
evolves. A very simple form for the differential equation can be given in this case by using 
(^), together with a harmonic trap potential Vt(x) = m(co'^x^ + Uyy"^ + ujyZ'^)/2; and the 

Thomas- Fermi approximation fijy = (^IbNucUxtOyCUzm^^'^ /16n\^"^ ^ , (with, however, a linear 
interpolation as ^ to give the correct noninteracting value of /i(0) = h{uJx + ujy + uJz)/2) 
yields a differential equation which can be easily integrated numerically. 

We present solutions for the parameters of the original rubidium ^ and sodium p| 
experiments in Fig.2. In both cases there is a latency time during which the condensate 



is initiated by the spontaneous term (the last term inside the curly brackets in (|23|) , after 
which the stimulated term (the term proportional to n inside the curly brackets in (PB|)) 
takes over, causing a rapid growth until saturation sets in when the condensate chemical 
potential /i„ approaches the chemical potential fi of the bath. 

The timescales for the growth of the condensate are of the same order of magnitude 
as experimentally observed, although no measurements have been published. In comparing 



5 



with experiment, one should bear in mind that this treatment (i) neglects fluctuation effects 
and (ii) neglects the quasiparticle effects. Inclusion of these should not change the growth 
curve in the region where n is substantially larger than 1, but could possibly speed up the 
process by which the first 100 or so atoms enter the condensate. 

The solutions in Fig. 2 assume that /i and T, the chemical potential and temperature of 
the bath of noncondensed atoms, are constant. They nevertheless exhibit the fundamental 
nature of the process of condensation. For a treatment more appropriate to comparison with 
experiment one must couple the condensate growth equation (^) to appropriate time devel- 
opment equations for the process of evaporative cooling, such as those of |T5|-pT[. There are 
three principal timescales in the problem; the timescale of equilibration of the noncondensate 
"bath", which is very fast, the timescale of condensate growth, as given by solutions of (p3|), 
and the timescale of the evaporative cooling process, which is in practice considerably slower 
than both of the others. Under these conditions one would expect that a model in which 
the noncondensate "bath" distribution function is considered to be always thermalized for 
particles below the "cut" energy, which we shall call rjkT, and is zero above this "cut" en- 
ergy would be valid. Provided the cooling process is slow enough, we can use (|2T| ) with the 
resulting time- dependent T and jl (values appropriate to the truncated distribution), after 
modifying the evaporative cooling equations to take account of the transfer of particles and 
energy between condensate and noncondensate. 

One should also note that the Boltzmann function with a cutoff at the energy r/fcT, 
typically with rj = 5-7, is a distribution which is significantly out of thermal equilibrium; 
for 7y = 5 or 7 we find respectively 12.5% and 2.9% of the full Boltzmann distribution is 
above rjkT. However collisions which are of the correct kinematic configuration to permit 
one of the atoms to enter the condensate are a selection of the full Boltzmann distribution in 
which the density of states factor, proportional to E"^ for the harmonic oscillator, becomes 
approximately independent of E; for them only 0.67% and 0.091% have energy greater than 
rikT. Thus the truncated distribution behaves like a genuine thermal distribution for the 
collisions which can populate the condensate. 

Using this more complete model of condensate formation, we have simulated a number 
of different evaporative cooling paths that are consistent with the published descriptions of 
Bose condensate formation (e.g 0,3). Although differing in detail, the results are broadly in 
agreement with those shown in Fig 2. One of the features of the experimental process is that 
the cooling process is halted, and the system allowed to thermalize for a short period before 
the condensate is observed. Our simulations show that the nonequilibrium noncondensate 
distribution evolves to a true Boltzmann distribution in a few mean collision times, and it 
is easy to show that the new temperature is always less than T, while the new chemical 
potential may be larger or smaller than /2. Typically, in the regime appropriate to the 
experiments, fi increases during this relaxation, and may change from a value below to a 
value above fic ■ This crucial step in the formation of the condensate, of fi evolving to 
exceed /ic, may thus have occurred only during the relaxation process. 

The condensate growth equation (^) is like the kind of equation one finds for laser. Thus 
there is the spontaneous emission term (the +1 inside the curly brackets), and a gain term. 
The gain here is determined entirely by the difference between the condensate chemical 
potential which is quantum mechanically determined, and that of the noncondensate 
"bath" /I, which is determined by statistical mechanics. The fact that the quasiparticles 
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play no significant role in the process is analogous to the behavior of a multimode laser, in 
which nearly all photons go into the mode with the highest gain, even if it is only marginally 
the highest gain. 

One can also conclude that adaptations of the condensate growth equation for different 
configurations of the noncondensate "bath" will prove a useful tool in the eventual design 
of an atom laser, or "Boscr" . 

This work was supported by the Marsdcn Fund under contract number PVT-603, and 
by Osterreichische Fonds zur Forderung der wissenschaftlichen Forschung. 
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FIGURES 



FIG. 1. Fig.l. The condensate and noncondensate bands 



FIG. 2. Fig.2. Condensate growth for (a) Rubidium and (b) Sodium. Scattering lengths are 
5.71 nm and 2.75 nm, respectively. 
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